Robust quantum enhanced phase estimation in a multimode interferometer 
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By exploiting the correlation properties of ultracold atoms in a multi-mode interferometer, we show how 
quantum enhanced measurement precision can be achieved with strong robustness to particle loss. While the 
potential for enhanced measurement precision is limited for even moderate loss in two-mode schemes, multi- 
mode schemes can be more robust. A ring interferometer for sensing rotational motion with non-interacting 
fermionic atoms can realize an uncertainty scaling of i/{Ny'7j) for particles with a fraction ri remaining 
after loss, which undercuts the shot noise limit of two mode interferometers. A second scheme with strongly- 
interacting bosons achieves a comparable measurement precision and improved readout. 
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Quantum interferometry aims at improving measurement 
precision with limited resources by way of quantum entan- 
glement of the constituent particles, say photons or ultra-cold 
atoms fU-lll]. However, entanglement is rapidly degraded by 
decoherence, e.g. due to particle loss. This makes it difficult 
to employ entanglement practically. Here we show how this 
limitation can be circumvented with robust multi-mode quan- 
tum superposition states. 

Interferometers estimate a parameter, e.g. a phase angle cf), 
by observing an interference pattern. A useful tool for quanti- 
fying the attainable quality of precision with a limited number 
of quanta N is the quantum Fisher information (QFI), Fq. 
This is the quantum analogue of the classical Fisher informa- 
tion and it is related to the uncertainty of the phase measure- 
ment via the Cramer-Rao bound, Scj) > I/^/Fq [4]. For 
unentangled particles (or independent single-quanta measure- 
ments) in the presence of loss, Fq = Nrj, where 77 is the frac- 
tion of particles remaining. This yields the shot-noise limit 
S(j) > l/^/Nfj. 

In two-mode interferometers, quantum-entangled input 
states can improve this to a maximum QFI of iV^, called the 
Heisenberg limit fUBl-lTll- This provides the lowest uncer- 
tainty allowed by quantum mechanics for linear interferom- 
etry schemes, i.e. when is a single-particle observable |8]. 
The Heisenberg limit is achieved with maximally entangled 
states known as NOON states lUlHHl- However, these are 
difficult to make |1^ and fragile to particle loss, reducing the 
QFI to Fq = N'^rj^ ifiol] . In the presence of only a small 
amount of loss (7/ < 1) the NOON state provides less pre- 
cision for large particle numbers than unentangled particles. 
Improving the robustness of entangled states is the subject of 
much recent work 1 10 - 16]. However, even when optimizing 
the input state for a given loss and particle number |lOl4l2r . 
the best achievable scaling is Fq ^ N for large N lIlTll . Thus 
Fq ^ N"^ is unattainable for two-mode schemes in the pres- 
ence of loss. 

This limitation does not apply to multi-mode schemes f\^. 
If sensitivity scales with mode number and arbitrary modes 
are accessible, there are no strict limits on the QFI for either 
classical (unentangled) or quantum input states. Unentangled 



states still suffer from a shot noise limit with the familiar scal- 
ing ^ N, when the particle number is varied in a given set of 
modes. Increasing the QFI can either be achieved by entan- 
glement within the given set of modes or by accessing more 
sensitive modes, while the latter procedure is the more robust 
against losses. The key questions in this context is how to pre- 
pare input states where mode occupation scales with particle 
number. For fermions it is natural to occupy many modes due 
to Pauli correlations, which forbid the occupation of a single 
mode by more than one particle. For bosons, interactions lead 
to similar coiTelations and the occupation of many modes. 

Here we demonstrate that a ring interferometer with non- 
interacting fermionic atoms can realize a QFI of Fq = N'^r/, 
which retains sub-shot noise scaling of the uncertainty in the 
presence of particle loss. Furthermore, we show that using 
strongly interacting bosons produces a similar precision with 
improved readout resolution. Quantum correlations are es- 
sential to enhance the sensitivity beyond the shot noise limit 
of an unentangled (classical) input state by providing access 
to a larger mode space rather than maximizing entanglement. 

Ring interferometer - Specifically, let us consider a ring 
interferometer with ultra-cold atoms to measure rotation an- 



gles ||19|]. Individual atoms, with an angular momentum hk, 
sense rotation by acquiring a phase shift fcc/). Here (/> = ijjtm 
is the rotation angle accumulated during the time interval tm 
and for a rotation rate lo. A single-particle interferometry 
scheme could estimate through comparing, by interference, 
the phase shifts accumulated by different angular momentum 
modes, say ki and ^2- Since the phase shift difference is 
(^2 — ki)4), we obtain a QFI of Fq = (^2 — ki)Nri in the 
case of Nrj successfully repeated measurements. By adjusting 
the angular momentum difference {k2 — ki)h, the measure- 
ment precision can be amplified without invoking entangle- 
ment, similar to Ref. [20]. We are going to show below how 
the natural properties of correlated multi-particle systems al- 
low us to scale the amplification factor with the particle num- 
ber in order to achieve a scaling beyond the shot noise limit, 
where Fq cx N"^. In contrast to NOON states, this scaling is 
not affected by loss. 

We now consider a ring interferometer containing N par- 
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tides. For the purpose of interferometry, the system may be 
prepared in a binary superposition 1^) = (l^i) + 1^2)), 

where \K) ~ X^n^^''^"!^) muhi- 
index n = {...,n^i,no,ni, ...) with fixed particle number. 



where aj, creates a particle with angular momentum hk. The 

notation Y^'"'^^ implies the additional constraint J2k ""^^k = K 
fixing the total angular momentum. The precise composition 
of these states depends on the interactions and quantum statis- 
tics of the particles involved, and has important consequences 
for the robustness properties. 

The superposition state is sensitive to rotation via the inter- 
action term Hn — —ojL, with the total angular momentum 
L = y^i. h katdk- This makes it suitable for rotation sens- 
ing il5ll2lll . After the time interval tm it evolves into 
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where an overall phase was ignored. 

We can quantify the ability of the quantum state to precisely 
estimate the angle </> by the QFI. It is independent of the mea- 
surement procedure and is given by 



(2) 



for a pure state, where |^''(</>)) = dl'^ {(f>)) / d(j) U2 

For the state given by Eq. ([Til we find Fq = {Ki — K2Y, 
which gives 6(t) > l/\Ki — K2\- To reach a scaling of 
Fq ^ N"^ we, therefore, require \Ki — = N. For a two 
mode system, where the modes differ by only one unit of an- 
gular momentum, this is only possible with the NOON state. 
This restraint does not apply for states described by more than 
two modes, because \Ki — = N can be achieved with 
many different configurations. As will be shown, this allows 
huge improvements to the robustness of sub-shot noise limited 
measurements. 

Particle Loss - For a system coupled to a zero temperature 
environment the evolution of the system in the presence of 
particle loss is described by the master equation |2JJ 
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where p is the density matrix, p — dp/dt, dk is the annihila- 
tion operator of mode k, and F is the loss rate, which is taken 
to be equal for all modes. It was shown in Ref. flT] that it 
does not matter whether the loss occurs before, during or after 
the phase shift is acquired. 

Equation ([3j is solved by describing p by iV + 1 density 
matrices, p^'^^^\ having the particle number N - 



N 



(N-){t)p(^-''^ 
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(^-'^+i)alistimeinde- 



Here - i V°' 

nere p - J^_^_^_^ 2^^, 

pendent and normalized to Trp'^""^' = 1, and g'^^^"^ (t) is 



a time dependent coefficient. Because the pt^"") operate on 
distinct orthogonal subspaces, the total QFI for a given loss 
rate is given by 



N 



(5) 



where F^ 
lows that 



is the QFI of [11]. From Eq. © it fol- 



(6) 



which has the solution g'-^ (t) 



-rt 



'{l — T])'^, where 

So Fq is easily determined for all 77 if the iV + 1 
values of Fq^'^ are known. We now compare the precision of 
unentangled atoms and NOON states with the precision of a 
fermionic and strongly-interacting bosonic superposition state 
in the presence of loss. 

Fermionic System - Consider an even number N — 2n of 
spin polarized fermions in the states | — n) = !!"=-« K^'-) 
and \n) = YVi=-n+i '^Ik^'')- Each one is a (shifted) Fermi 
sea with total angular momentum —nh and nh, respectively. 
The superposition \tpF) = ") + !"-)) thus has a QFI 

of N"^ and realizes sub-shot noise limited scaling. The state 
\^pF) can also be written as 
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71 — 1. 



(7) 



■i— — n+l 



which shows that, effectively, only a single particle partici- 
pates in the superposition. While the states | — n) and \n) are 
degenerate ground states of the kinetic energy / {2mE?) 
at zero rotation, the superposition \iPf) emerges as the only 
ground state when the degeneracy is lifted in a weak external 
potential that breaks the rotational symmetry. 

To evaluate the performance under loss, we consider the 
density matrix after removal of v particles 



P 



(TV!) 



ri-l 

E 

fcl<...<fc„-l = -n+l 
1 n-1 



afc„-i-afci|x)(x|aii-«L_i 



fci<...<fe„ = -n+l 



(8) 



where \x) = nr^T-n+i '^I 1^^'-) ^ filled Fermi sea with N— 1 
particles. 

Since all terms are represented on mutually orthogonal sub- 
spaces, their QFI can be calculated individually and summed. 
We note that the first summation has no off-diagonal terms, 
so the QFI is simply zero. The second summation is a sum of 
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FIG. 1. Top: A schematic showing the implementation of the scheme 
in time. For a description of each step see the main text. Bottom: A 
visual representation of the system. Ultra-cold atoms are confined to 
a ID optical ring potential with a rotating barrier. The measurement 
precision is determined in the presence of particle loss. 



(^^ ^) pure states, which all have QFI of N'^. This leads to 



p(JV-i/) _ N-i> 
— N 

for a given loss, 



N^. Substituting into Eq.|5]gives the QFI 



Fq - N^Tj. 



(9) 



Significantly, the QFI is iV^ for no loss and decreases at the 
same rate as unentangled atoms when loss is considered. This 
is the expected result, because the state is a single-particle su- 
perposition with a difference in angular momentum of Nh. 
The role of the remaining atoms is to facilitate the creation of 
this superposition in the ground state. 

Although this scheme gives an excellent precision and ro- 
bustness to loss, the readout would require distinguishing the 
difference in momentum of a single particle, which is diffi- 
cult. Furthermore, an even number of particles in the initial 
state needs to be ensured (e.g. by post-selection on the basis 
of total angular momentum) in order to exploit its degeneracy. 
We now consider a scheme using strongly interacting bosonic 
atoms with improved readout options that is also insensitive 
to the even-odd particle number parity. 

Tonks-Girardeau System - Bosonic atoms in a tightly con- 
fined wave guide at low densities and strong repulsive inter- 
actions realize the Tonks-Girardeau (TG) gas of impenetrable 
bosons 1 24]. The TG gas maps one-to-one to non-interacting 
spin polarized fermions |25l], and thus has similar properties. 
The non-rotating ground state, however, is non-degenerate 
and thus the previous scheme needs to be modified. A binary 
superposition with angular momentum difference of Nh can 
be created with the help of a rotating narrow barrier potential 
as discussed in detail in Ref. |26]. Here, we calculate the QFI 
and suggest an interferometry scheme that can take advantage 
of the theoretical precision bound. 

Specifically, we consider bosonic atoms in a one- 



dimensional ring trap of circumference L = 27ri? at zero tem- 
perature. They are stirred by a barrier of strength 6, which 
rotates with angular velocity lu — hVt/ {mL"^), where J7 is a 
phase induced around the ring in the co-rotating frame (see 
Fig ID- The system is described by the Hamiltonian 
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ki ,k2,q—^oo 



(10) 



where Eq = 2tt^ fi^ / [mL'^) is the smallest non-zero kinetic 
energy of a single atom and g > is the interatomic interac- 
tion strength. 

With external rotation at = tt and 6 = 0, the ground state 
is degenerate between the states |0) with zero and |A^) with 
Nh angular momentum, regardless of the interaction strength. 
A finite barrier strength < 6 <C g^/N /2 lifts the degeneracy 
and It/)^) = (|0) ± |iV))/-\/2 become eigenstates, where l-i/)^) 
is the ground state. This state has a QFI of N"^ and thus real- 
izes sub-shot noise limited scaling. In the TG regime where 
formally g oo, the Bose-Fermi mapping |25J allows us to 
map \tp^) onto a fermionic state similar to Eq. 

A measurement scheme for rotation is shown schematically 
in Fig.[T] Step 1 corresponds to the creation of the initial state, 
where uj is slowly increased up to wq = nh/ {mL'^), corre- 
sponding to f2 = TT (for details see Ref. ll26ll ). Once the initial 
state has been created, the state of the system will not change 
until the barrier's rotation rate is altered as indicated by step 
2. The additional rotation to be measured, Aoj, is then non- 
adiabatically applied (step 3), inducing an additional phase 
of Ail around the ring. The Hamiltonian describing the sys- 
tem is then H' = EZ-oo Eo{k-^ + 4^)' where 
we ignore the last two terms in Eq. ( fTOl i. because the barrier 
height is small and interactions do not couple states of differ- 
ent total angular momentum or change | K) when the rotation 
is changed. The system is then allowed to evolve for time tm 
as shown by step 3. This establishes a phase difference be- 
tween the two parts of the superposition and thus mixes \ipg) 
and iV'ij)- At this point we calculate the QFI. 

A possible readout scheme is to remove Aw non- 
adiabatically leaving the barrier rotating at its original rate, uq 
(see step 4 in figure[T]|. The rotation of the barrier is then adi- 
abatically reduced to a point where the states \tjjg) and \tjjg) 
evolve into |0) and |iV), respectively (step 5). The trapping 
potential is then removed and the atoms are imaged. The |0) 
state will have a peak at the center of the image while \N) 
will have a dip due to the different momentum distributions. 
By adiabatically reducing the interaction strength 5 to a small 
value before removing the confining potential, the states |0) 
and \N) can be transformed into (ao)^|vac) and (a|)^|vac), 
respectively. This allows for an efficient distinction of the two 
outcomes. The QFI can be shown to saturate the classical 
Fisher information for the total angular momentum. 

The ring interferometer scheme can, with appropriate mod- 



4 



Tl 


1 1 

\ 1 — Femiions 




— ^ — Tonks-Girardeau 




- — — Optimal 2-mode ■ 

\ NOON 

V — — Unentangled 











0.2 0.4 0.6 0.8 1 

FIG. 2. The uncertainty of (f) — Aujt^ for different fractions of parti- 
cles remaining, rj, and N — 5. The lines show 5<l) for non-interacting 
fermions and TG (M = 18) superpositions compared with the opti- 
mal 2-mode (as given in 1 10]), NOON and unentangled states. The 
non-interacting fermions and TG (AI — 18) superpositions afford 
better precision than the unentangled particles and the optimum two- 
mode initial state for all loss rates. 



ifications, also measure rotational phases with a NOON state 
[(^o)^+('^i)^] I vac), which is obtained for small interactions, 
or with a state of unentangled particles [(aj + a|)^] |vac), ob- 
tained for b > gVN /2. 

Using Eq. (|5]l we now compare the effects of particle loss 
(during time <,„) on the precision capabilities of the three 
states. The spread out momentum distribution of the TG 
state is different from fermions and limits numerical sim- 
ulations of Eq. ( [Tol l currently to 5 particles, where rescal- 
ing of g ensures accurate results in a truncated basis of 18 
momentum modes 1261 IztIi . The QFI for a mixed state is 
Fq = TT[p{(j))A'^], where p{(f)) is the density matrix of the 
system and A is the symmetric logarithmic negativity, defined 
as = ^[Ap{(t)) + p{(j))A] m. Intheeigenbasisof p(0) 
this is {A)ij — 2[p' {(f>)]ij / {Xi + Xj), where Xij are the eigen- 
values of and p'{(f>) — dp{(t>) /&((). If + Xj =0 then 

The results are shown in Fig. |2]for = 5. As expected 
the precision of the NOON, fermionic and TG superposition 
states are equivalent when there is no loss, ?7 = 1, and the pre- 
cision of the unentangled state is much worse. As the loss rate 
increases, however, the unentangled state soon outperforms 
the NOON state. Importantly, the fermionic and TG superpo- 
sition states outperform the unentangled state for all loss rates 
and therefore could prove extremely valuable for metrology. 

For comparison with previous work, we also show the op- 
timized initial two-mode state of Refs. ifiol [l3ll for N — b 
in Fig. 121 We see that whilst it outperforms the NOON state, 
for stronger loss its precision gradually approaches that of an 
unentangled state as the fraction of particles remaining, 77, de- 
creases. Importantly the precision of the fermionic and TG 
superposition states are better than the precision of the opti- 
mised two-mode initial state for all loss rates. Not only this. 



the optimal two-mode initial state very much depends on the 
amount of loss as the structure of the state changes with 77 
thereby making its experimental implementation difficult. 

Conclusion - We have shown that robust sub-shot noise lim- 
ited measurements are made possible by strong correlations 
of fermionic atoms due to the Pauli exclusion principle, and 
in TG systems due to strong interactions. The proposed states 
not only offer improved scaling of measurement precision, but 
also outperform optimized two mode states for small particle 
numbers. This is of significance for metrology as any state of- 
fering an increase in precision over unentangled states has the 
potential to alter the way precision measurements are made. 
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Massey University for hospitality. 
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